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CALCULATION OF SLOW ENERGY NEUTRON SPECTRA 

GJMarchukf GoAJljasova, VM Morozov, V,V ^Smelov, VnAModakov, 

There are three most /essential aspects in the-neutron thermalization problem: the develop- 
ment of the theory of neutron-matter interaction and the evalution of corresponding constants, 
needed for calculation; the development of the computational algorithms for solving of 
mathematical problems with needed accuracy; and, at last, the problem of the interpretation of 
slow neutron flux spectral data. On the whole, an the present paper the questions of second 
aspect of the'problem will be considered* 

The-slow energy neutron scattering problem was formulated by Hurwitz and Cohen [1,2]. 

The -most complete solution of this problem was obtained in the theory of neutron scattering by 
nuclei of monoatomic gases [1—4]. In the last years the-necessity turned out of a more 
careful analysis of neutron scattering mechanism, taking into account the molecular binding 
effect. The -results of these -investigations were reported in detail at the Brookhaven 
Conference on Neutron Thermalization [5]. 

It should be noted then that in the last years some important experimental data on the 
neutron thermalization problem have been obtained [5-7]. After these investigations the-more 
careful comparison of experimental data with the theory based on different physical models is 
possible* It should be -also observed that progress of our knowledge of physical processes in 
nuclear reactors makes the development of more perfect and exact methods of the solution of 
neutron transport equation evident. 

The -peculiarity of the-problern of calculating the thermal neutron distribution lies in the 
fact, that the low energy neutrons not only lose their energy, -but acquire it. Due to this 
effect, the integral operator in the transport equations is Fredholm-type operator. The 
transport equation may be solved by spherical harmonics method in Pji-approximation, 
by Sp-inelhod, by the method of characteristics, by Monte Carlo method, and others [8—15]. 

The survey of some mathematical methods of solution of the transport equation, -their 
confrontation and comparison of the theory with the -experiment [6—7] are -given in present 
paper. A spe^'ial attention is paid to the-problern of calculating neutron flux angular distribu- 
tions cello 
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I. CALCULATION OI' NEUTRON FLUX AND NEUTOON IMPORTANCE BY THE 
SPHERICAL HARMONICS METHOD. 

The stationary slow neutron flux in the case of a cylindrical cell may bc'described by the 
following integro-differential transport equation: 

sin(9(cos^ M -SiM. ^)+a(r, \)4,(T,v,e,ib) - 
or r Sift 

n \ 

-/ dv’/dfl’cis(r;v’-iv, €l’4n)0(r, v',d’,((r')-q(r, v) 

where agfr; v’-« v, n ’40) is the differential scattering cross-section; a(r, v)«as(r, v) + aa{r, v) 
is the'tptal cross-section; d and ifr are the meridional and azimuthal angles, respectively (Fig.l). 
Let the function <^(r, V, bean the form 

( r . V, -p(v)^(r ; v, (2) 

where p(v) is chosen to take a more complete. account of the neutrpn flux (^(r, v, 6, 
character of variation in variable v. The energy spectrum of formally homogenized cell was 
used as the . function p(v). From above one may conclude) that v(r, v, 6, if/) is slightly 
dependent on the. velocity v over the entire. thermalization range ■0<v<Vgroup. 

Divide the range 0<v<vgroup into m intervals vj<v<Vj+i (j-1,2 ... , m). Within each energy 
group Vj<v<vj^.i the.function vtr.v, S,il/) can be approximately regarded as independent of 
neutron velocity. Then, substititting (2) in equation (1) and integrating the latter over 
everyone of the. groups, we. obtain the. following system of equations: 

sin^(cosv>^ - 2iS± )+«()) (r)^(i)(r, 0,^)_ 

dr f dif/ 

m /n ... 

where ^(j>(r. d,,^r)-/i+l^(r,v,d)dv , q{j)(r)-/j-^lq(r,v)dv . 

' vj 

a(j)(r)-i f a(v)p(v)dv , p;-/'^ p(v)dv 

* Vj 

as'*'(r, fi'4n).,;L / ®‘^^p(v’)av’ as(r, v’4v, fl’-«i)dv, 

pc e I 


j • 

Let us expand 'the < function Qg in Legendre 'polynomials 


«s ’(r,lT’-^)-i. 2^^ 2|±L aJ">(r)P„(go) 
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ond neutron flux into spherical harmonics, so that 
e;, ^ 1 ,( 1 ) (r)Pj^(cos 0 )+ 

nil Ll Sr Wcosnu^ Pi™) (cosfl) 

P 3 : is the ' approximation, when all (r) of n>4 are neglected* Because of the. scattering 
function (r,ft'-n) independence of neither medium absorptivity nor of its size* in series 
(3) one may take 'into account only finite number of terms without connecting this fact with the 
order of approximation. For practice. in (3) only two terms can be. taken into account. As 
a result for each energy group we. obtain tlie system of six differential equations for (|)|^ (r), 
which in a matrix. form is^) 


a — +pTor + 2 (i)=So 
ai.^+LTiq,+ 2 ir. 
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S^i](r)-4,rq(i){r)+(l-Sjg) 1 $(|) (r)ai"i(r) 
Sg(r)-<l-aji)|^ <!,(•) (r)aH 


®Of ® 1 * 'Toi Tl, So, are. some matrix coefficients [10, 16]. 
Excluding the .vector r from system (5), we. obtain the. equation: 

+ ,^k 4J£ -(3alSo + l-AKl,- 

dr 2 ^ r 2 

=-3viSo^-3y(^ + LSi)-iWSi 

In equation ( 6 ): Sq, ^ are diagonal metrixes, 

5a, 5a/6] ■•[0,0, 4h] 


(5) 


( 6 ) 


x) 


For the sake of simplicity, an index of the.group number is omitted. 
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I=j.-r-^ — , = •» +i — - 

(a I g, Z 1 a 

Equations (6) are to be completed by the boundary conditions, providing symmetry of 
neutron flux at the'Cell centre-and reflection condition at its boundary [10, 16]. ) 

The solution of sys^m (6) is obtained by Zeidel iteration method, In this method the. 
functions and SjY at n»th iteration step are calculated in the following schemer 

S®n - + J.. ‘I>oo!n-l«ol-*i al-*i 

The. solution of the. multi-group theory must satisfy the.neutron balance relation: 

5 Rgroup ^y, | (r) rd, (7) 

)-l 0 J.l 0 

If the. iterating is not accomplished, balance relation (7) is not achieved. If one. trying to 
obtain balance relation (7) during the. iteration process, it is possible. to hasten considerably 
the. convergence. of the. process [17]. In present paperitwas made by the.sourcerenormalization. 

Each:onegroup equation of system (6) is solved by the method of the matrix. factorization. 

The. functions (-)-/,;S(i)(r, 0,(j)dr. thus obtaine^a, ,r’y be used for the. calculation cf 
effective. thermal constants, needed for further calculations of critical reactor parameters, in 
some. cases, however, it is important to. know the angular dependence. of the neutron-flux 
^(i)(r>^,^^). in some. experimental investigations of thermal neutron spectra tl.e neutron beams 
of the. certain directions are usedi Formula (4) reveals the . possibility to calculate. the 
neutron fluxes in any given direction. Below we cite the expressions of neuttpn fluxes with 
different velocity directions (see.Pig.l). 


(r, C (r) + 54>2ii* (r) ] , 

0(i){r, 0r. ^ ^ « 1_ (r) - 5 (i>2i} (d 

<|i(i)(r, „ ]) =1_ [<l>oo (r)-|. ( 1 ) 2^0 W + : 

± (34.iY(r) -.7 +7^ q,^^^r))]. 


After the-neutron flux <;6(r, v, 0, \jj) the • neutron imp>ortance*^*(r, v, in relation to some 

process is of interest. Forexamplei one can say about the < neutron importance -in relation to the 
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neutron absorption in block, etc.[18-20]. The-neutron importance is a solution of the adjoint 
equation, 'which in multigroup representation has the -form: 


dr r dij/ 

Jl / 0 *(l)(r, ig' (r, UM') do* -q*(j)(r, 0 ,^) 

1=1 


(9) 


The concrete 'form of the riglit side of equation (9) is connected- with physical process in 
relation to which the neutron importance is determined. In particular, for the calculation of the- 
neutron importance in relation to the block capture* one-must put: 


iV . I 4'^ w in the block 

beyond the block 


The solution cf system (9) can be obtained by the same methods, as for system (10. 

The neutron importance in the 'combination with the neutron flux distribution function is 
widely used in the 'perturbation theory [18*21]. 

11. CALCULATION OF NP:U1'?10N FLUX BY .^^lETHOD. 

To calculate'the thermal neutron flux in a reactor cell, the-Sn-method, proposed by 
Carlson [12,13], may be used successfully. In this section we -shall discuss the -application of 
this method to multigroup system (D, the -anisotropy of scattering will be taken into account 
only in transport approximation. 

The solution of the corresponding system of equations with the -reflection conditions at the- 
cell boundary is realized-by the-following three types of iterating: inelastic transition 
iteration, elastic scattering iteration and boundary condition iteration. The-equation iterated 
is of the form: 


yl^yZcos^ <^</>y|tvP ^^s,*t-,p 

dr ^ 


dt/f 

where j is the energy group number (j*l, 2 ....,m), y«cosd, 


(j) 




-J? fi'-'lriNjijw+qOjCr) , 


rr 0 

the -total transport cross-section is 

= /3i-»kr) + aj^(r) 

tbc-SGUttcring cross-section is 


(10) 


( 11 ) 

( 12 ) 

(13) 

(14) 


(r) . 
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The.functions ajf (r) and (r) are. step functions, which are some constants for each 

velocity group; q()l(r) is the linear function in the intervals, limited by the knot points of the 
variabler. 

The Iterating of the inelestic tmnsitions of S-number consists in the successive solution of 
equations /lO/ for all j with the.functions Q^'^r), being known from the ptevious iteration 
Step of inelastic transitions. 

The Iteration, of the elastic collisions of t-number consists in the. solution of the system 
/lO/ for given .j with the.function (r). obtained in the.previous iteration of the elastic 
scattering. The.function (r) as a result of all iterations of elastic collisions is found. 

And, at last, .by iterating the boundary conditions the.function 03^1 which satisfies the 

reflection boundary conditions: V/ v) ^ ^ • i*. • j rvu 

^ IS obtained. There is a 

condition at the-outer boundary of a cell in the. iterating of the boundary conditions: 

(IS, 

in all 7r/2<^<rr interval. 

T'he. computational formulas of the Sn-method, needed for solution of the iterated 
equations (10), were. obtained by the general for all geometries scheme, reported in [U]. The 
approximate. solutions of the multigroup system obtained by these formulas satisfy the 
neutron balance condition. In accordance. with this scheme.the range of variation of 
variables r. y is covered by a- net of knot points yp).The.knot points ^ the 

interva [0.1] are chosen as coinciding to Gauss integral formula; tlie.knot points in the 
mterval [0,,r] are.chosen as equidistant with the.w/n step (the number n coincides the 
n-approximation); the knot points rk in the.interval [0,Rgroup] are chosen arbitrarily 
bu tb.condition should be observed, that the points of discontinuity of the functions at?(r), 

(r) and.q(])(r) were among - points. 

xhe integral operators and ^^rdr are. successively applied to the equations of 

system (10) at y=yp with the assumption of linearity ^(i) in ^ and r in the intervals 

.q->ipnp,e..b,ai„, „bich co».ec. ,he. 

Addiog lo this spslsm the.e,u.tip«s, obtained after sohstitnting in (10) «.i,h the 

sobseqng.. oppiioauo. of the oper.«4^_^r ,e hove the.eo*^^^ 

diacnssod .bore req.ies , of iterations. 

To reduce a nnmbet of iterations, »e.used the comspondios to„,n|as of the acceleration 
of convergence. of iterating. acceleration 

The.convergence.of boundary ^pnditions iterating was hastened by extrapolation of the 

of iterations r (Li.st.mil. ee.bti "“»■>- 

Ti.e.con.ergence.of elastic collision iteralions .as hastened 1, the es.lntatloa method of 
...^r..,«.de.,..,„„s 114,251 Theddeaef this ...hod is i.the appro^ate d.^i” “of the 
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iterated function deviation from the exact solution of the equation by it^ variation after one 
Iteration step. In our case the corresponding formulas can be obtained in the following way. 
The solution of equation ( 10 ) <^3 (r,(fr,y) can be represented' in the form: 

where ' A 0 s^(,(r,(i-,y) is the. iterative. deviation, corresponding to the iteration of t-number. 
Integrating (17) over 0 - and y, we get 

(r)- (r) - ANg^iJ, (r) ( 13 ) 

vjrhere * ^ ^ 

ANs5/(r)=L / r A,^i;J,(r.dr,y)dy n 9) 

" 0 b 

Using (17) and (18), we. obtain from (10): 

Vi^ cos(ir _Yl-y2 sinji. dAd-i^t +aHV) aJ') - 

=di-‘i(r)ANs^^\^l (r): 

from which the. equation for the iterative deviation Aijiglt, follows: 

- 3i-»i(r)ANj| +Ag 3 (j? (r) 
where (i) • • 

^'Is,t. -Ns,t.(r)] ( 22 ) 

The boundary conditions of this equation are the previous reflection ones. 

The next step is the • estimation of the iterative deviation value. ANg'l (r) or^ in other 
words, the. approximate solution of equation (21). From the fact, that with -jig 
the function A(;is|t becomes more. smooth, we. obtain the equation for an estimation of 

iterative -deviation A 0 g^^yr,^,y): 

A/IZ;2cos#i^ £^ 5 ^ +.!!’wA4i -4qS W; ( 23 ) 


No iterations of the elastic collisions are needed. for the. solution of this equation. 

The. approximate. function of average flux* thus obtained: 

^sjt W = Ngjf (r) - ANg^(,(r) ( 24 ) 

ANs^^\^r)-L;d^/A0i;J, (r,^fr,y)dy „ (25) 

IS used then as a new approximation in iterated equation (10). All iterating and introducing of 
iter;=tive corrections last till the ratio ANg (! (r) for all r becomes less than the. value <. 

S A ^ * 


) 
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As the function gives only approximate value of the iterative deviation from 

the 'exact solution <{>s\ then to reduce numerical work, it must be determined in more low-order 
approximations in comparison to the approximations, used in calculation of itself* 

In cases when the conditions for introducing of the correction from (24) were not held* the 
convergence of elastic collision iterating was hastened by the factorization of the iterated 
function by the'balance normalization factor, The same acceleratipn of convergence was 
applied to inelastic transition iterating. 

III. MONTE CARLO CALCULATION OF NEUTRON FLUX. 

Using the-probability characteristics of neutron Interactions with nuclei and those of the 
cell geometry, the neutron wandering in the cell is modelized by means of a digital computer 

machine* At the cell boundary a neutron is reflected -mirroriy with no velocity change* The 
source-in the thermalization range is constructed with the assumption, that for velocities 
v>Vgroup a neutron is scattered (elastically and isotropically in the centre mass system) by 
fixed free nuclei and that at v>Vgj.oup the neutron spectrum is Fermi spectrum. 

All probability characteristics of neutron interactions can be tubulated, but it needs a large 
volume of an operative computer memory. 

In the case -of a gas model for moderator the -neutron scattering by nuclei of a certain type 
can be modelized as follows. The-absolute value of the-velocity Vq of nuclei with which 
neutron collides, as selected. Fot this purpose -in relation to the following equation is solved 

^2^^4()')+3'A2(y) , 

< 

“0-^4 — <;i 3 (/S)] +/3[<^i(y) — y>l3 

K= iiz - 

^2 ^4(y)+302(^)'+|t2’“'^3(/3)]+/3[|--<Al(/3)3 


where ^- V/V-r, r-Vo/VT. V is the 'neutron velocity, -Vx is the most pro babiye nuclear 

velocity at the -temperature T^K, K is a random number, 0<K<1, te^"*' dt. The 'functions 

** 0 

0ji(z) are to be computed beforehand* 

Then the -angle ^ between the-neutron motion and the nucleus as a target direction is 
selected: 


cos^«^ [K(l+a)3+(l_K) ll_a|3] f. 

a=*Vo/V. The-second-angular parameter is determined from the condition of equivalence of 
azimuthal directions. From the known neutron and nuclear velocities the velocity of a 

scattered neutron is calculated-by 


[V+MVo+Mlf-^ole 

M+1 


e is a random isotropic unit vector. 
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It is known a number of ways to store the information on the flux during the neutron wander- 
ing [25]. As a rule, the spectrum in the moderator is rapidly extablished. The difficulties are 
connected with greatly absorbing domains of a small volume^ The calculation of the flux in 
any part of the cell may be significantly hastened, if one takes into account the contribution of 
each collisiony after which a possible path of neutrcn scattered crosses its volume* In the 
cells with moderator volume predominance the computations can be significantly accelerated 
when replacing the moderator regions, remote from the block, by the corresponding surface 
sources. The characteristics of sources (spectrum, angular distribution eta) are to be 
computed beforehand* in onedimentional cells it is realized in an especially simple way. 

The graphite-uranium cell was computed by this method (page 13). 

As a rule, to get the exactness of computation, practically accepted, (7-.10).10'^ collisions 
are sufficient. The results of computations, in which more than lO^ collisions were selected, 
are shown in Figs 12,13. As for spectrum in graphite (Fig.l3b) further corrections were made* 
Below the approximate values of standard errors of some -characteristics of the fluxes 
computed are represented (in %): 


Cell 

Zone 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Uranium-water 

Block 

0.89 

0.83 

2,2 

1.4 

2.4 

5.9 

3.6 


Moderator 

0.14 

0.36 

0.94 

0.55 

1.1 

2.3 

2.1 

Uranium-graphite 

Block 

0.73 

0.89 

2.9 

1.3 

1.9 

2.4 

3.5 


(1) - total flux, (2) — total absorption, (3), (4), (5), (6), (7), - fluxes of variable x=: V/ ^21 ^ q 
in intervals (OJ), (1.2), (2.3), (34'), (4. 8), respectively, To=r300°K- 

In computations three-order generator of pseudorandom numbers of a known type: 

(mod M) was used. We -put p=55, M-224, X^^l and Xo=3. The generator was tested by the 
computation of Maxwellian spectrum in a cell, with no absorption taken into account, and by 
Kolmogorov criterion -cuS at the different sample -sizes. 

The results are satisfactory. 

It is useful to smooth the density functions, obtained by statistical method. The following 
way is very useful too [26]. The empirical function of probability distrtbution F(x), cor- 
responding to the -calculated density function, is considered. At the knot point Xj^, 4c«s, 
s+1, , n-s (n+1 is the-total number of knots) F(x}^) is replaced by P(xk), where Ps(x) is a 

polynomial of certain degree, which minimizes the value S [F(xi^^i -P(xk+i)]^- By 

differentiating a more smooth distribution function, thus constructed, the smoothed density 
function is obtained. We-used the polynomialsof the fifth degree at Sa8/n«80/. 

IV. RESULTS OF CALCULATIONS 

The-described above methods were used for the calculation of the -space-energy slow neutron 
distribution along the cells of uranium-water and- uranium-graphite lattices, studied 
experimentally by Mostovoi [6,7]. 
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The uranium-water assembly is a set of blocks of natural uranium 3.5cm in diameter, located 
in the sites of a trianpjular lattice witli the 5 . 5 cm spacing. 

The uranium-graphite assembly consists of the same uranium blocks, located in the sites of 
a square lattice with 20 cm spacing. In the experiments the uranium blocks are separated by 
an aluminium layer from the moderator. 

In calculations the experimental hexagonal and square cells were replaced by the equivalent 
circular threezone cells. 

Now we proceed to describing the calculations of cells in P 3 - approximation. The 
calculation of slow neutron flux was performed in the energy range 0<K<^]g^Q^p = 0.67 ev in 
15-group approximation. The sources of thermal neutrons caused by slowing down were 
calculated with the assumption, that at E>Eg^.Q^p the neutron spectrum is Fermi spectrum and 
molecular binding effect at E>EgJ.Q^p is not of importance. The inelastic molecular scattering 
and neutron interaction with the -crystal lattice at 0 <E<E^Q^p were taken into account on the 
basis of Turchin results [22,23]. 

For fissionable isotopes the deviation of the capture cross-section from ^/v low was 
taken into account. 

The space distributions of average thermal neutron flux are given in Figs 2,3. 

dy/<^(r,v, 0 ,^)df 2 , 

and also directed fluxes /see ( 8 )/; 

‘I'llW-/ = 0)dV, ci)|°) (r)«/^(r,v,^..E, V'-0)dV, 

o X o 

^ ^ ^ V ^ ^ V 

(r)-/ \ 

It is seen from Figs 2 and 3, that neutron flux <l>ii(r) (neutron velocities are parallel to a cell 
axis) is ^'eaten out’ * much stronger in the block, than the average flux <5>oo(r). From physical 
point of view, this results from the fact, that the capture probability of neutrons moving in 
a block along its axis, is greater, than for neutrons of other directions. As the distance from 
block increases, becomes larger, than ^I>oo(r), and this fact can be also simply 

explained from physical point of view. 

Let us discuss in detail curve -3 (Figs 2 and 3). It describesa neutron beam, directed along 
a cell diameter. The left branch of this curve (r<0) corresponds to the neutrons, moving 
towards the cell axis, » and right one does to the neutrons moving in an opposite direction 

<^^)(r). In the -most part of its path in a block, the neutron beam loses, as a result of capture, 
more neutrons, than acquires from other beams due-to scattering. This results in a monotone 
character of curve 3 . in a block. The disturbance of the beam monotony when leaving the 
block is connected with the general increase of neutron density near the block boundary. 
Crossing the-aluminium layer, the -beam <l^®)(r) does not considerably change due to absorption 
and scattering, and the abrupt decrease of curve 3 in a layer is due to its geometry. 
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The distribution along diameter of those . neutrons, velocities of which are. perpendicular to 
an axial cross-section of the cell, <J)^’'/2)(r) is represented by curve. 4 in [•ig.2. It is evident 
that at r»0 <!^°)(0)=.(I^'r/2)(0)B<l^n')(0). * 

Now we proceed to an analysis of the thermal neutron energy spectra in cells /Figs 4-9/. 

The abscissa is plotted in variable 1/x (x-v/\f^ .T-dOOOK). The ordinata is the neutron 
density n{l/v)-v2n(v)-v,^(v). In Figs 4,6,8 the neutron, spectra at the centre of the cell are shown 
The curve 1 represents the. neutron energy distribution, averaged over the angles 

j-/vgi(r=O,v,0)dn. The curve 2 corresponds to neutron beam energy distribution along the 
cell axis v.^(r-O,v,0»O) and curve. 3. presents the spectrum of the. neutron beam perpendicular to 
the. cell axis v,4(r-0,v,^-|.). From physical considerations spectrum 2 is to bcharder than 
global spectrum 1, but spectrum 3 - softer. 

It is well in agreement with the calculations. The. Spectra given are compared with 
experimentally measured ones. The experimental method was such that the. beams extracted 
from the blocks were mainly of neutrons, directed along the. axis of the cell. It is natural 
therefore to compare the experimental spectra with the fu.nctions v,^(r-0,v,^-0). 

However, one should expect that the experiraental spectrum will be a bit softer than the 
theoretical one> since the. neutron beam extracted has the neutrons of other directions. In the 
uranium-water cell this prediction is confirmed by the calculation. In the. uranium-water cell 
t IS effect is weak due to a considerable softeningof the spectrum calculated in the moderator 
as compared with the experiment. The origin of this softening is not yet clear. 

The comparison of the.same experimental spectrum with the global spectrum 1 fv<^(r4),v,(?,,A)dn 
reveals that the experimental one.is harser. In Figs 5,7,9 the comparison of the tlfe'^oretical and^ 
exporunental spectra at the cell boundary is given. Calculations show the neutron flux at the 

ceh boundary IS practically isotropic. FigslO, 11 present space and energy distribution of the 

neutron importance. m relation to the absorption in a block for the. uranium-water cell The 
presenceof maxima in the.curves 1 and 2 (Fig.ll) in therange 3.0<x<4.0 is caused by ‘the 
resonance in the. absorption cross-section ,u235. n is of interest to confront all the three 
methods given above* For this purpose. the- thermal neutron fluxes in uramum-water and 
uranium-graphite.cells were calculated by each method. The. computations were. made on the- 
basts of gas^model, stnce.in this case the. scattering indicatrix can be selected by Monte 
Carlo method completely. The results of the.computations ate.shown in Figs 12 and 13 The 
mutual attachment of the spectra calculated by different methods was made in the block as 

well as in the.moderator at the maximum. There is a good agreement between the spectrl given 
within the errors of the methods. ® 

We.ar=,™.ch obliged Dr. VJ.lloseovoi .„d Dr. V.S.Dikarev for .heir i„.ere,. i. the-pre.r.. 

work and useful discussions during its reaKzation. 
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